Localization of coherent exciton transport in phase space 
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We study numerically the dynamics of excitons on discrete rings in the presence of static disorder Based 
on continuous-time quantum walks we compute the time evolution of the Wigner function (WF) both for pure 
diagonal (site) disorder, as well as for diagonal and off-diagonal (site and transfer) disorder In both cases, large 
disorder leads to localization and destroys the characteristic phase space patterns of the WF found in the absence 
of disorder 
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I. INTRODUCTION 

Ever since the emergence and development of quantum me- 
chanics, there has been a major interest in the crossover from 
quantum mechanical transport to the corresponding classical 
transport. However, there are far-reaching differences in the 
usual mathematical description of the two different processes. 
In classical physics, the phase space is spanned by conjugate 
variables, such as position and momentum, whose time de- 
velopment leads to classical transport. Quantum mechanical 
processes, on the contrary, take place in Hilbert space. One 
approach to overcome these differences was already presented 
in the early days of quantum mechanics, more than 70 years 
ago, by Wigner BjJ, |2[|. He introduced a function, now known 
as Wigner function (WF), which is a (quasi) probability in 
a quantum mechanical phase space spanned by position and 
momentum variables. The WF is a real valued function and 
in this respect compares well with the classical Boltzmann 
probability distribution in phase space. However, it is not al- 
ways positive. WFs and related phase space functions, like the 
Husimi function, are widely used in Quantum Optics 0,0] but 
also for describing electronic transport, see e.g. HIE01- 

Apart from the quantum-classical crossover there are also 
extremely interesting purely quantum mechanical phenom- 
ena. For example, Anderson has shown that there is no quan- 
tum diffusion for some random lattices (gl d, OH, an ef- 
fect nowadays called (strong) localization, where the quan- 
tum mechanical transport through the lattice is, in essence, 
prohibited by the potential energy surface. Anderson's hop- 
ping model for electron transport has also turned out to be 
useful in describing excitons in disordered systems, e.g., 
El EES QUI. 

In fact, Anderson's model is closely related to the so-called 
continuous-time quantum walks (CTQWs) with disorder. Re- 
cently, it has been shown that the motion on a graph, described 
by CTQWs, can be exponentially suppressed by the disorder 
flol . Originally, (unperturbed) CTQWs were introduced in 
the context of quantum information as the quantum mechan- 
ical analog of continuous-time random walks (CTRWs) H 1 711 - 
Here, the underlying, discrete connectivity of the structure on 
which the transport takes place determines the Hamiltonian. 
Since CTQWs also model exciton transport over various dis- 
crete structures 11811 . they are closely connected to other ap- 
praches to study (coherent) exciton transport phenomena on 
discrete graphs, for instance, in the contexts of polymer 11911 , 



atomic [20] or solid-state [21] physics. 

In this paper we consider the dynamics of excitons on dis- 
crete rings under static disorder, focussing on a quantum me- 
chanical phase space approach. We compute the correspond- 
ing WFs numerically and study the effects of both pure diag- 
onal (site) disorder, as well as of diagonal and off-diagonal 
(site and transfer) disorder on the (coherent) transport. The 
WFs are then compared to the unperturbed case, for which an 
analytical treatment is possible B22h . 

The paper is organized as follows. In Sec. m we briefly re- 
view how we model coherent exciton transport on graphs and 
define the appropriate Hamiltonian for rings with two types of 
disorder. After introducing the discrete WF in Sec. Jill we lay 
down the procedure of our calculations in Sec.|lVl The subse- 
quent Sees. IVl to I VIII show the numerical results for coherent 
exciton transport on rings with disorder. We close with our 
conclusions. 



II. COHERENT EXCITON DYNAMICS 

In the absence of disorder, the (coherent) dynamics of exci- 
tons on a graph of connected nodes is modelled by the CTQW. 
Here, the CTQW is obtained by identifying the Hamiltonian 
of the system with the (classical) transfer matrix, H = — T, 
see e.g. 11711 (we will set h = 1 in the following). The transfer 
matrix of the walk, T = (T/j), can be related to the con- 
nectivity matrix A of the graph by T = —7 A, where for 
simplicity we assume the transmission rates 7 of all bonds to 
be equal. The matrix A has as non-diagonal elements Aij the 
values — 1 if nodes I and j of the graph are connected by a 
bond and otherwise. The diagonal elements Ajj of A equal 
the number of bonds fj which exit from node j. 

Now, the states \j) associated with excitations localized at 
the nodes j span the whole accessible Hilbert space to be con- 
sidered here. In general, the time evolution of a state \j) 
starting at time to = is given by \j;t) = exp(— iHt)\j). 
By denoting the eigenstates of H by \$>q) and the eigen- 
values by Eq the transition probability reads 7rij(t) = 
Kl|exp(-im)|j)| 2 = \^ e eM-iEeW\*e)($e\j)\ 2 . 
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A. Dynamics on rings without disorder 

The unperturbed Hamiltonian for such a CTQW on a finite 
one-dimensional network of length N with periodic bound- 
ary conditions (PBC) takes on the very simple form H°|j) = 
2\j) — \j — 1) — \j + 1), where we have taken the trans- 
mission rate to be 7 = 1 and j = 0, 1, . . . , N — 1. The 
eigenstates of the time independent Schrodinger equa- 
tion H |$$) = Eq\^ q ) are Bloch states, which can be ex- 
pressed as linear combinations of the states \j), see Refs. 
|22i|23|] for details. 



B. Dynamics on rings with disorder 

Now we introduce (static) disorder by adding to the unper- 
turbed Hamiltonian H° a disorder operator A, i.e., by setting 
H = H° + A. We let the disorder matrix A = (Aij) have 
non-zero entries only at the positions for which Hij 7^ 0. 
For different strengths of disorder, the elements Aij = Ajj 
are chosen randomly (drawn from a normal distribution with 
the mean value zero and the variance one), which we multiply 
by a factor of A, which then can take values from the interval 
[0,1/2]. Note that under these assumptions for the (static) dis- 
order the connectivity of the graph is essentially unchanged, 
i.e., there are no new connections created nor are existing con- 
nections destroyed. Therefore, the only non-zero matrix ele- 
ments of H are those of the initial A. The action of the new 
Hamiltonian H on a state \j) reads thus 

H|j) = (H + A)|j) = 2|j)-|j-l)-|j + l) 

+ 2A jd \j) - Ajj-i\j - 1) - Ajj+i\j + 1). (1) 

In the following we consider two cases of disorder: 

(A) Diagonal disorder (DD), where Ajj ^ and Aij = 
for I / j. We assign a random number to each Ajj, a 
procedure which leads to N random numbers. 

(B) Diagonal and off-diagonal disorder (DOD), where we 
choose a random number for each Ajj and for each Ajj-i. 
Thus, 2N random numbers are needed. 

Introducing disorder into the system in this way has conse- 
quences for the relation between the CTQWs and the CTRWs. 
In CTRWs the transition rates, given by the entries of the 
the transfer matrix T, are correlated, i.e., for each site the 
sum of the non-diagonal rates for transmission from it and 
the diagonal rate of leaving it are the same. In the cases 
considered here, a direct identification of the Hamiltonian H 
with a classical transfer matrix T is, in general, not possi- 
ble anymore. However, the DOD and DD Hamiltonians are 
widely used in quantum mechanical nearest-neighbor hop- 
ping models, to which also the CTQWs belong. Furthermore, 
we still consider transport processes on graphs which have 
the connectivity matrix A, but the direct connection between 
H and T is lost. Of course, one can maintain this connec- 
tion by imposing constraints on the A, e.g., by requiring that 
2Ajj = Ajj -1 + Ajj+i for all j. 



III. WIGNER FUNCTIONS 

The WF is a quasi-probability (in the sense that it can be- 
come negative) in the quantum mechanical phase space. For 
a phase space spanned by the continuous variables X and K, 
the WF reads 111 HH 

W(X,K]t) = ^ J dY e iKY (X -Y/2\p(t)\X + Y/2), 

(2) 

where p(t) is the density operator. In the following we only 
consider pure states, i.e., p(t) = | j; t) (j; t\. Note that the WF 
is normalized to one when integrated over the whole phase 
space. 

For a discrete system having N sites on a ring, where we 
choose to enumerate the sites as 0, 1, . . . , N — 1, (x\j;t) is 
only defined for integer values of £ = 0,1,..., TV — 1, and the 
form of Eq.© has to be changed from an integral to a sum. 
Now the WF resembles a discrete Fourier transform, which 
requires N different k- values. These k- values may evidently 
be chosen as k = 2ttk/N, again having ft = 0, l,...,iV — 1. 
According to I22I1 we use, for integer x and y, the following 
discrete WF 

1 N ~ 1 

Wj{x,K-t) = - elkV (x-y\j;t)(j;t\x + y). (3) 

y=0 

The marginal distributions are now given by summing over 
lines in phase space, e.g., we get 

N-1 N-1 N-1 

N-1 

= ^2 s y,o(x-y\j;t)(j;t\x + y) 

y=0 

= \(x\j;t)\ 2 =7T xJ (t). (4) 

We note here that there are forms similar to E g. Q] ) in defining 
discrete WFs, see, e.g., Refs. 

In order to compare to the classical long time behavior, we 
define a long time average of the WF by 

T 

Wj(x,K) = lim ^ / dtWj(x,K,;i). (5) 

T^oo T J 


As for the marginal distributions for the WF we ob- 
tain marginal distributions for Wj(x,K) upon integration 
along lines in (x, n) -space. Especially, summing over k, 
J2 K Wj(x, k) = Xxj, one obtains the long time average of 
the transition probability n X j (t). 

IV. WFS ON RINGS 

A. WFs on rings without disorder 

The unperturbed Hamiltonian for a CTQW on a ring is 
given by H°, see Sec. Ill Al In this case, the discrete WF can 
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be calculated analytically by using a Bloch ansatz. It is given 
by HI 



N-l 



Wj(x, «; t) = ex P [ - i2 < 2n + K )( x ~ fil N ] 



n=0 



x exp { - z2£[cos(2tt(« + n)/N) - cos(27ra/iV)]}. (6) 

As shown in Ref. 12211 . the marginal distributions of 

Wj (x,K,;t) and Wj (x, k) are correctly recovered. 



B. WFs on rings with disorder 

Under disorder the Bloch property is lost, then the previous 
analytic approach does not apply anymore. We hence com- 
pute numerically the WFs for the corresponding quantum dy- 
namics of excitons on rings with disorder by using the stan- 
dard software package MATLAB. The differences in the fol- 
lowing computations are only due to the different Hamiltoni- 
ans H, depending on the specific type of disorder. 

Having the eigenstates \ $q) and the eigenvalues Eq, the WF 
is obtained by expanding Eq. © as 



N-l 



y=0 



x (x - y\$e>)($e>\j)(j\$e)($e\x + y). (7) 



Figure \T\ shows the WF according to Eq.© for N = 101 
and two different realizations of H with A = 1 /2 at time 
t = 40. Clearly, the details of the WF differ much for dif- 
ferent realizations. However, a trend of the effect of disorder 
on the dynamics is already visible: The disorder prevents the 
excitation to travel freely through the graph. Instead, there is 
a localized area about the initial site x = j = 50. We will 
study this in much more detail in the remainder. 
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FIG. 1 : (Color online) Snapshots of two realizations of the WF for 
N = 101 and DOD with A = 1/2 at time t 40. 

There appear non- vanishing values of the WF at nodes op- 
posite to the initial node j, i.e., at x w j + N/2, even at 
infinitesimal small times t. These values are not related to 
disorder but are due to the PBCs. This can be inferred from 
Eq. ([3]). Consider first the case of even N: for t = and x = j 
the only contributions to the sum in Eq. ^ are those for which 
y = 0. However, for t = and also for x = j + N/2 the WF 
does not vanish: There are non-zero contributions to the sum 
fory = N/2, since we have (j\j)(j\j + N) = = 1. 



Non-zero values of the WF for x = j + N/2 continue to 
show up at later times, see also [22] . The argumentation for 
odd N is similar: Here, however, non-vanishing values at 
x w j + N/2 start to appear as soon as the wave function 
(x\j; t) spreads over more than one node. The difference be- 
tween these non-vanishing WF- values for even and for odd 
N at short times tends to zero as 1/N for N large. We re- 
mark that the appearence of these non-vanishing WF- values 
for short times at x ~ j + N/2 does not lead to a finite 
transition probability 7Tj + jv/2,j ;(£), since at small t the sum 
in Eq.© practically vanishes. This is in line with the situ- 
ation for open boundaries, where the transition probabilities 
7r j+N/2,j(t) a l so vanish at short times. There, however, also 
the WF- values at x w j + N/ 2 are practically zero at short 
times. A thorough study of the influence of different bound- 
ary conditions on the WF is beyond the scope of this paper 
and will be given elsewhere. 

Once we have the WF, we also compute its long time aver- 
age. Now it is preferrable not to first evaluate Eq. © and to 
perform then the computationally expensive time integrals of 
Eq. ©, but to proceed at first analytically, obtaining 

w^k) = ^Y, eikv Y.~ 5 ^- E ^ 

y e,e' 

x {x-y\§ e >){®e>\3){3\®e){$>e\x + y), (8) 

with 8{E e - E e >) = 1 if E e = E Q > and S(E e - E e >) = 
otherwise. Eq. {3]) is in general much more accurate and 
computationally cheaper. 



C. Ensemble averages 

In order to have a global picture of the effect of disorder 
on the dynamics, we will consider ensemble averages of the 
WFs. For this we calculate the WF for different realizations 
of H and average over all realizations, i.e., for R realizations 
we compute 

1 R 

(WifaKrfU = -J2 [Wj(x : ^t)] r , (9) 

r=l 

where [Wj (x, t)] is the WF of the rth realization of H. 



V. THE ROLE OF DISORDER 

For all cases of disorder, we study the WF for different 
strength of the disorder, i.e., for different A with A ranging 
from A = 1/40 to A = 1/2. We exemplify our results for 
odd (N = 101) and even (N = 100) numbered rings, where 
in all cases the initial condition (t = 0) is a state localized at 
node j (in our two examples j = 50). Note that the WF of a 
localized state at x = j is equipartitioned in the ^-direction, 
i.e.,^(j>;0) = l/7V. 
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A. Diagonal disorder (DD) 

We start by considering diagonal disorder for a graph with 
N = 101. Figure [2 shows snapshots of (Wj(x,K]t))R for 
different values of A at different times. 
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FIG. 2: (Color online) DD: Ensemble average of WFs of the quantum 
dynamics on a ring of length N = 101 for A = 1/40, 1/10, 1/4, 
and 1/2 [columns (l)-(4)], each at times t = 1, 10, 20, 40, 100, and 
500 [rows (a)-(f)]. The initial node is always j = 50 and the average 
is over R = 1000 realizations. Red regions denote positive values 
of the averaged WFs, blue regions negative values and white regions 
values close to 0. The colormaps are always chosen to be the same 
for each row but might differ in different rows. The maximal val- 
ues of (Wj(x, K;t))R are denoted by small black regions; these are 
highlighted and exemplified by the arrows in panel (4f). 

The first column shows (Wj (x , t)) r for A = 1/40 at 
times t = 1,10,20,40,100, and 500 [Fig. |2la)-(lf)]. For 
this quite weak disorder, the patterns in phase space are sim- 
ilar to the unperturbed case, where "waves" in phase space 
emanate from the initial site x = j = 50 and start to interfere 
after having reached the opposite site of the ring (see Fig. 3 
of i22h ). However, at longer times differences become visi- 
ble, Fig O If). The pattern for the unperturbed case is quite 
irregular but with alternating positive and negative regions of 
the WF of approximately the same magnitude. Here, the dis- 
order causes a decrease of (Wj(x, k; t))n for x close to the 
initial site j = 50 and k- values in the middle of the interval 
[0, TV — 1] compared to the values of n close to or N — 1. 

Increasing the disorder parameter A, the patterns change 
profoundly. The wave structure gets suppressed and for all k 
a localized region forms about the initial site j = 50 already 
for small disorder (A = 1/10) and short times (t = 20), see 



Fig. [22c). 

For even larger values of A, for all k the formation of a lo- 
calized region about j = 50 becomes even more pronounced. 
Already for A = 1/4 this localized region forms for times as 
short as t = 10, see Fig. [23b). At A = 1/2, (Wj(x, n;t)) R 
stays localized for all times [Fig.[24a)-(4f)]. Also here, values 
of the WF at about k « N/ 2 are decreased, whereas values of 
the WF for k about the interval borders and N — 1 remain 
rather large for x ~ j, as indicated by the thin black region 
(see also the arrows), e.g., in Fig.[24f)- We further note that 
at high disorder the localized averaged WF is always positive, 
i.e., all fluctuations, present for small disorder, have vanished. 
We recall that the WF is normalized to one when integrated 
over the whole phase space. 
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FIG. 3: (Color online) DD: Same as Fig.[2]but for N = 100. 



Having an even number of nodes in the graph does not al- 
ter the picture drastically. By comparing Fig. [2 to Fig. [3j 
which shows the quantum dynamics on a ring of N = 100 
nodes with DD, one sees that the corresponding panels are 
rather similar; the localization effect of the disorder on the 
system stays the same. However, there are also differences 
at long times, compare (le-f) and (2e-f) in each figure. An- 
other difference between even and odd N lies in the fact that 
for even N there remains a region of alternating values of 
(Wj(x, K\t))n at about x = j + N/2, even for large disor- 
der, see column (4a)-(4f). This is due to the periodic bound- 
ary condition and, therefore, to the equal number of steps in 
both directions starting from any site on the network needed to 
reach the opposite site of the network. Obviously, the disorder 
does not destroy this symmetry. 
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B. Diagonal and off-diagonal disorder (DOD) 
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FIG. 4: (Color online) DOD: Same values of t, A, and R as used in 
Fig.0 

The second type of disorder is the diagonal and off-diagonal 
one (DOD), where both the diagonal and the off-diagonal ele- 
ments of the unperturbed Hamiltonian are randomly changed, 
as explained in Sec. Ill B I Figure @] shows (Wj(x, ft; for 
the same values of t, A, and R as in Fig. [2] Here and as for 
DD, the symmetry with respect to the line in phase space at 
x = j = 50 remains intact. Comparing to the DD case, there 
is a slightly weaker suppression of the waves in phase space. 
This effect is rather small but can be seen by comparing, e.g., 
the values of the WF in Fig.[2{2f) to the ones in Fig.[4{2f), note 
the different corresponding colorbars. However, the overall 
localization effect is very similar for both types of disorder. 
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evaluated in this way for TV = 101 at time t = 500 and for 
A = 1/40, 1/10, 1/4 and 1 /2. Also here, the localization ef- 
fect is clearly seen. Nevertheless, the details of the WF differ 
from those obtained for the other types of disorder. 



VI. MARGINAL DISTRIBUTIONS 

Integrating the WF along lines in phase space gives the 
marginal distributions. Since we saw that the effect of local- 
ization does not depend on the particular choice of the type of 
disorder, we display the marginal distributions only for DOD. 
Figure[6] shows the marginal distributions for N = 101 at time 
t = 100 for different A. Obviously, the details of the WF are 
lost when integrating along either the x or the k direction. 
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FIG. 6: (Color online) DOD: Marginal distributions of 
(Wj(x, ft; t)) R , for N = 101 at t = 100, for A = 1/40,1/10,1/4, 
and 1/2. (a) summing along the k direction, (b) summing along the 
x direction. 

While the transition probabilities, obtained by summing 
over all k [FigHa)], clearly show the effect of localization, 
this is not the case for the marginal distribution obtained by 
summing over all x [FigHb)]. In the latter case, one can- 
not readily distinguish the situations corresponding to differ- 
ent values of A. 



VII. LONG TIME AVERAGES 



FIG. 5: (Color online) Constrained diagonal and off-diagonal disor- 
der, see text for details: Averaged WF for N = 101 and t = 500 
[cases (f) in Figs.[2]to|4l, but for the same values of A and R as used 
in Fig. [2 



Another possibility is to constrain the disorder by having 
+ Ajj+i, which maintains the connection of 
H to the classical transfer matrix T. Figure \5\ shows the WF 



In has been shown in 12211 that for unperturbed coherent 
exciton transport and odd N (superscript °), most points in 
the quantum mechanical phase space have a weight of 1 /N 2 , 
namely we have 



W°Ax, k) = < 



l/N 2 ft^Oandanyx 
1/N k = and x = j 
else. 



(10) 
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For even N, the limiting WF reads 

{2/N 2 k 7^ 0, k even and any x 
l/N k = Oandx = jJ + N/2 (11) 
else. 

We note that Eq. <Qj} differs from Eq. (19) of Ref. JH, which 
is not correct. The long time averages of the WFs for even N 
are somewhat peculiar, since values different from zero appear 
only for even k, whereas the WFs themselves have values dif- 
ferent from zero at arbitrary times for all k. A numerical check 
(which we do not include here) for a finite line of N nodes 
without disorder, shows that these stripes in the long time av- 
erage in the present study are due to the periodic boundaries. 
For even N there are constructive interference patterns in the 
transition probabilities, since the number of steps in both di- 
rections is the same, see also Ref. l23ll . For a finite line with 
even TV, there are no stripes in the long time average. Of 
course, for short times and close to the initial node, the WFs 
for the line and the circle coincide, since the boundaries have 
no effect on the initial propagation. Nevertheless, one also 
has to bear in mind that the long time average is not a real 
equilibrium distribution. 

As a technical note, we remark that changing the order of 
the time and ensemble averages can lead to a considerable 
speed-up of the numerical computation. We checked numer- 
ically that the time average and the ensemble average indeed 
interchange, i.e., we have 

T 

{Wj(x,K)) R = (^lim^/ dtWj(x,K;t)) R 


T 

= lim^ J dt(Wj(x,K;t)) R . (12) 
o 

Now, using Eq. ® reduces the computational effort further. 
For computational reasons, we further interchanged the sum- 
mation over y with the ensemble average, i.e., 

{Wj( X , k)) r = (^E e%ky E ^ - E <>') 

y e,e' 
x(x- y\&o')(®e'\j)(j\®o)(®o\x + y)) 

v e,e> 
x{x-y\^){^\j){j\^ e ){^e\x + y)) . (13) 

Again, we carefully checked that Eq. (IT3l) yields the same re- 
sult as Eq. (fT2l) . 

Now, the disorder changes also the limiting WF quite dras- 
tically. Starting from high disorder of A = 1/2, we expect 
from Figs.[2];4a)-(4f) to Il4a)-(4f) that the long time average 
of the averaged WF will look basically the same. Figure [7] 
shows the limiting averaged WF for N = 101 and DOD ac- 
cording to Eq. (fT3l) . Indeed, for large A, the limiting averaged 
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FIG. 7: (Color online) DOD: Limiting averaged WF, (Wj(x, k))r, 
for N = 101 and A = 1/40, 1/10, 1/4, and 1/2 [panels (1) to (4), 
respectively], according to Eq.([T3l). 



WF is comparable to the corresponding averaged WF, com- 
pare Figs.[7t4) and0t4f). Close to the initial node x = j = 50 
and there along the ^-direction, (Wj(x, k))r has large (posi- 
tive) values for k about and N — 1 which decrease by going 
toward k = N/2. Here, the minimum of this decrease de- 
pends on the particular type of disorder. 

Also for small A there are significant differences to the 
case without disorder. From Fig. Hi) we see that the dis- 
order "smears out" the localized value Wj(j,0) = l/N 
in the case without disorder, see Eq.JTOl). Specifically, the 
(x = k = 0) -value decreases while the neighboring ones 
increase. For A = 1/10, the onset of localization about 
x = j = 50 is already seen, see Fig.[7t2), and becomes more 
and more pronounced as A increases, Fig. [33). Furthermore, 
all other values of (Wj(x, k))r for x ^ j decrease with in- 
creasing disorder, as can be seen by the decreasing size of the 
light pink region, which corresponds to values close to 1 /N 2 . 
In fact, the values of the limiting averaged WF for x distant 
from x = j = 50 drop to zero, shown as white regions in 
Fig. [7] 
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FIG. 8: (Color online) Same as Fig.[7]but for N = 100. 
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For even N we found that without disorder the limiting WF 
shows a peculiar "striped" distribution, caused by the PBCs, 
see Eq.dTTl). By switching on the disorder, these peculiarities 
vanish. Figure [8] shows the limiting averaged WF for DOD 
and N = 100. Although for small A there are some remain- 
ders of stripes left, these disappear completely for higher val- 
ues of A, compare Figs.[8tl)-(4). Note further that the second 
peak of Wj (x, k) at x = j + N/ 2, see Eq.fTTlL transforms for 
increasing disorder to an oscillatory line in the ^-direction at 
x = j + N/2. 
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FIG. 9: (Color online) DOD: Marginal distribution J2 K (Wj (x,k))r 
for different A and (a) N = 100 and (b) N = 101. 

By summing now again along the ^-direction, this pecu- 
liar behavior vanishes and also for even TV we get a localized 
marginal limiting probability distribution at x = j. Figure [9] 
shows the marginal distribution (Wj (x, k))r for even and 
odd TV with DOD. As expected, for even A", Fig. Oa), there 
are no remainders of the peculiar distribution of (Wj (x,k))r 
at x = j + N/2 for large A. Furthermore, for high disor- 



der, the marginal distributions of (Wj(x, k))r have a shape 
similar to that of (Wj (x, t))R, compare Figs.Oa) andOb). 

In general, the difference in the long time average 
J] K (lfj(i, k))r between even and odd A" strongly depends 
on the disorder. Without disorder, the two cases are distinct 
ll22ll . For odd N the constructive interference at the node 
j leads to only one peak in ^Z K (W 3 -,(x, k))r, see Eq. (20) 
in l22n . However, for even N the constructive interference 
at nodes j and also at node j + N/ 2 leads to two peaks in 
y }2 K/ (Wj(x J k))r, see Eq. (21) in Ii22l1. With increasing dis- 
order, the second peak in ^2 K (W 'j(x, k)) r for even A" at 
x = j + N/2 vanishes and gives rise to only one peak at 
x = j, see Fig.[9j Nonetheless, for large but finite Af the WF 
itself [given in Eq.©] still allows us to distinguish between 
the odd and the even case, see Figs. [7] and 
VIII. CONCLUSIONS 

We have analyzed the effect of static disorder on the coher- 
ent exciton transport by means of discrete Wigner functions. 
We have studied numerically the dynamics on a ring of A^ sites 
in the presence of pure diagonal disorder and also of diagonal 
and off-diagonal disorder. The previously found characteris- 
tic patterns of the unperturbed WF in the quantum mechanical 
phase space are destroyed by the disorder. Instead, the WF 
shows strong localization about the initial node. Integrating 
out the details of the time evolution by considering the long 
time average of the WF, shows an even more pronounced lo- 
calization. 
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